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SN 2009ib: A Type II-P Supernova with an Unusually 
Long Plateau 
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ABSTRACT 

We present optical and near-infrared photometry and spectroscopy of SN 2009ib, 
a Type II-P supernova in NGC 1559. This object has moderate brightness, similar 
to those of the intermediate-luminosity SNe 2008in and 2009N. Its plateau phase is 
unusually long, lasting for about 130 days after explosion. The spectra are similar 
to those of the subluminous SN 2002gd, with moderate expansion velocities. We es¬ 
timate the ®®Ni mass produced as 0.046 ± 0.015 Mq. We determine the distance to 
SN 2009ib using both the expanding photosphere method (EPM) and the standard 
candle method. We also apply EPM to SN 1986L, a type II-P SN that exploded in the 
same galaxy. Combining the results of different methods, we conclude the distance to 
NGC 1559 as D = 19.8 ± 3.0 Mpc. We examine archival, pre-explosion images of the 
field taken with the Hubble Space Telescope, and find a faint source at the position of 
the SN, which has a yellow colour {{V — I)o = 0.85 mag). Assuming it is a single star, 
we estimate its initial mass as Mzams = 20 Mq. We also examine the possibility, that 
instead of the yellow source the progenitor of SN 2009ib is a red supergiant star too 
faint to be detected. In this case we estimate the upper limit for the initial zero-age 
main sequence mass of the progenitor to be ^ 14 — 17 Mq. In addition, we infer the 
physical properties of the progenitor at the explosion via hydrodynamical modelling of 
the observables, and estimate the total energy as ^ 0.55 x 10®^ erg, the pre-explosion 
radius as 400 Rq, and the ejected envelope mass as ^ 15 Mq, which implies that 
the mass of the progenitor before explosion was ^ 16.5 — 17 Mq. 

Key words: supernovae: general - supernovae: individual: SN 2009ib - galaxies: 
individual: NGC 1559 
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1 INTRODUCTION 

Type ITP supernovae (SNe II-P) are by far the most com¬ 
mon among the core-collapse SNe representing about 60% 
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of them iILi et al volume-limited sample). They can be 

classified by the strong hydrogen lines in their spectra and 
by their characteristic light curves that show a near constant 
Inminosity during the first ~ 100 days of evolution. 

SNe II-P with lower absolute magnitudes, slower 
expansion velocities and smaller ®®Ni masses than the 
majority have been studied extensiv ely. SN 1997D was 
the first sub l umino us SN analysed (iTuratto et al.l 1 19981 : 
|Benettiet_alJ [2001 ), and was followed by e.g. 1999br 
l!Hamuvll200ll : Pastorello et al.l l 2004h. 2002gd (ISpiro et al] 


|2014). 2003Z iutrobin et al.ll20()7l:ISDiro et al.ll201.^ r2005cs 


I Thkats_&_^inkL^^ 

I Fraser et al J|2Qlll) 


Pastorello et al. 2' 


and 2009md 


Recently, a few analyses of SNe ITP that have parame¬ 
ters falling in between those of the normal and subluminou s 
objects have been pu blished. SNe 2008i n (IRov et ab 201lh, 
2009js llGandhi et al.ll^l3h . and 2009N llTakats et al. 12014 1 
showed mixed properties, similarities to both groups. The 
observed properties of SNe 2008in and 2009N, such as their 
luminosities, spectra, expansion velocities were almost iden¬ 
tical. By estimating the progenitor para meters at the explo¬ 
sion for SN 2009N, iTakhts et al.l (l2014| j found that the ex¬ 
plosion energy and the estimated ®®Ni mass of this object are 
between the typical values of the subluminous and normal 
SNe II-P. iGandhi et ^ (l2013l l showed that the characteris¬ 
tics of SN 2009js place this object in between the sublumi¬ 
nous SNe and SN 2008in. The studies of these intermediate 
SNe suggest that subluminous and normal SNe ITP do not 

form distinct groups. _ 

_ Similar conclusions ha v e bee n drawn bv lAnderson et ^ 

ll2014h and ISanders et al.l (l2015l l. who carried out studies 
that examined t he lig ht curves of a large number of SNe. 
I Anderson et al.l (12014 ) co l lected T-band light curves of 116 
SNe, while Sanders et al. Ii liolsl) used gri 2 j/-band photome¬ 
try of Pan-STARRSl objects. Both groups concluded that 
the light curve properties of SNe II, such as the absolute 
magnitude, plateau length, decline rate have a continuous 
distribution. 

SNe II-P are thought to emerge from stars with main se¬ 
quence masses of > 8 Mq . This is supported by cases, when 
red supergiant progenitors with masses between 8 — 16 Mq 
have been directly detected on archival, high- resolution pre¬ 
explosion images (recent examples include ISmar tt l2009l: 


[Fraser et al.ll20l'il : IVan Dvk et al.l[2012l : iMaund et al. 


2014ll . 


Despite these studies, however, we still do not fully under¬ 
stand the causes of the observed diversity of these objects. 
The analysis of new and different objects can help to better 
understand their progenitor and the explosion mechanism. 

In this paper we present data and analysis of an ob¬ 
ject, that seems to belong to this group of “intermediate” 
SNe ITP, but also shows significant differences from those 
previously studied, further blurring the line between the sub¬ 
groups of SNe ITP. SN 2009ib was discovered by the Ghilean 
Automatic Supernova Search ( GHASE) project on August 
6.30 UT 2009 in NGG 1559 (jPignata et al.l [2009ll. This 
galax y have hosted thr ee other SNe, 198 4J ( Evans et al.l 
I1984 t ype II), 1986L llEvans et al.l Il986l . type ITP) and 
2005df (IWang fc Baadell2005l . type laj. SN 2009ib was clas¬ 
sified as a type ITP SN similar to S N 2005cs a week after 
its explosion on 2009 August 9.1 UT (iStritzinger &: Folatellil 

This paper is organised as follows. In Sect. [2] we present 



Figure 1. The field of NGC 1559, the host galaxy of SN 2009ib. 
The stars used for photometric calibration are marked (see Tables 
lAll - IA^ . The image was taken on 2009 October 17 with the 
PROMPT5 telescope in V band. The insert shows the positions 
of all the SNe discovered in this galaxy, SN 1984J was a Type II, 
1986L a Type ITP, 2005df a Type la SN. 


optical and near-infrared photometric observations, estimate 
the extinction towards the SN and examine its colour evo¬ 
lution. In Sect. [3] the optical and NIR spectra are shown 
and compared to those of other SNe ITP. We estimate the 
distance to SN 2009ib in Sect. |4]using the expanding photo¬ 
sphere method (EPM) and the standardized candle method 
(SGM), and apply EPM to SN 1986L, another SN ITP that 
exploded in the same host galaxy. In this section we estimate 
the time of explosion of SN 2009ib as to = 55041.3 MJD 
(2009 July 29.3 UT), the epoch we use throughout the pa¬ 
per. We calculate the bolometric luminosity and the pro¬ 
duced ®®Ni mass in Sect. [S] In Sect. O we examine the 
pre-explosion images of the SN site taken with the Hubble 
Space Telescope (HST) searching for the progenitor, and 
we also determine the progenitor properties via hydrody- 
namical modelling of the SN observables and compare the 
obtained physical parameters with those of other SNe ITP. 
We compare the light curve and spectra of SN 2009ib to 
the non-LTE (non-local thermodyna mic equilibrium) time - 
dependent radiative-transfer models of iPessart et al.l (l2013ll . 
Finally, in Sect.[7]the main results are summarized. 


2 PHOTOMETRY 

Optical photometry was collected using multiple telescopes 
with UBVRI and u'g'r'i'z' filters, covering the phases be¬ 
tween 13 and 262 days after explosion. The basic reduc¬ 
tion steps of the images (such as bias-subtraction, overscan- 
correction, flat-fielding) were carried out using the stan- 
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dard IRAI0 tasks. The photometric measurement of the SN 
was performed using the point-spread fnnction (PSF) fitting 
technique via the SNOOPY packag43in iraf. The calibration 
of the photometry was carried out by observ i ng standard 
fields llLandoltlll992l : [Landolt fc Uomot3l2007l : ISmith et al.l 
l2002l i on photometric nights. With the help of these images, 
magnitudes for a local sequence of stars (Fig. [T] Tables HU 
and HU) on the SN field were determined and used to cali¬ 
brate the SN measurements (Tables (Adi and TASI) . 

Near-infrared photometry was obtained using the Rapid 
Eye Mount telescope (REM) in JH bands. Dithered im¬ 
ages of the SN field were taken in multiple sequences of 
five. The object images were dark- and flatfield-corrected, 
combined to create a sky images then the sky images were 
substracted from the object images. The images were then 
registered and combined. Photometry was carried out with 
PSF-fitting, also using the SNOOPY package. In order to cal¬ 
ibrate the measurements, magnitudes of local sequence of 
stars were obtained f rom the 2MASS catalogue (Table HU 
ISkrutskie et al.ll2006li . The JH magnitudes of the SN can 
be found in Table IA6I and in Fig. [2] 

SN 2009ib was discovered ~ 8 days after explosion 
(±2 days. Section 14.1. Ill , at the beginning of its plateau 
phase. On the plateau the magnitudes in Vr'Ri'Iz' bands 
are roughly constant. In u'Bg' bands the magnitudes de¬ 
crease more rapidly in the first ~ 40 days, while in the JH 
bands the brightness increases for ~ 75 days (Fig. [21. 

The plateau phase of this SN is unusually long, longer 
than that of most SNe II-P. After the plateau the brightness 
drops almost 2 magnitudes in ~ 30 days. Unfortunately, 
the tail phase is less well-sampled, with the last measure¬ 
ment at 262 days. This phase is powered by the radioac¬ 
tive decay of ®®Co to ®®Fe. The expected decline rate, as- 
suming full trappi ng of the gamma-rays, is 0.98 mag/lOOd 
llPatat et al.l Il994l l. Between days 160 and 262 we mea¬ 
sure the decline rate in the bolometric light curve (SectjSj 
as 1.12 ± 0.05 mag/lOOd, while in V band as 0.78 ± 0.05 
mag/lOOd, in R band as 1.14 ± 0.09 mag/lOOd. Note, that 
number of data points to fit the light curve tail is low in 
every band and our observations only cover the beginning 
of the tail phase, therefore these measurements are quite 
uncertain. 

Fig. [3] compares the shape of the light curves of a 
sample SNe II-P. We included in this sample the sublumi- 
nous SN 2005cs, the normal SN 2004et, and the interme¬ 
diate SN 2009N, which had a brightness similar to that of 
SN 2009ib. Some important parameters of the SNe used for 
comparison throughout the paper can be found in Table [T] 
Fig. [3] shows that the main differences in the light curves 
are the length of the plateau and the brightness of the tail 
relative to the plateau. SN 2009ib has a very long plateau 
and a relatively small drop to the tail. SN 2009N has a sig¬ 
nificantly shorter plateau, while the drop to the tail is about 
the same as for SN 2004et. The normal SN 2004et and the 


^ IRAF is distributed by the National Optical Astronomy Obser¬ 
vatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under the cooperative agreement 
with the National Science Foundation. 

^ SNOOPY is a package originally designed by F. Patat and later 
implemented in IRAF by E. Cappellaro. 



Figure 3. Comparison of the BVRl light curves of SNe II-P. 
The light curves were vertically shifted so their peak magnitudes 
match those of SN 2009ib in each band. The explosion epochs 
of the comparison SNe are very well-constrained (Table (TJ. The 
diversity of the shapes is reflected mostly in the length of the 
plateau and the drop to the tail. SN 2009N (red) has the same 
brightness as SN 2009ib (black circles), however its plateau is 
shorter. SNe 2005cs (green) and 2004et (blue) have about the 
same plateau length, but their absolute brightness differ by more 
than 2 magnitudes and the drop in the end of the plateau of 
SN 2005cs is greater. 

subluminous SN 2005cs has similar plateau lengths, while 
the drop to the tail is greater for SN 2005cs. 

In order to compare the plateau length of SN 2009ib to 
that of other SNe, we determined the epoch of the middle of 
the transition from the plateau to the tail ( fpr) via fitting 
the analytic function of lOlivares et al. l ||2010|) to the V band 
light curves of the SNe. Note, that [Anderson et al.l (l2014h 
used the same method on a much larger sample to estimate 
tpT- The difference between tpT and the explosion epoch 
(to) of the comparison SNe can be found in Table [U In the 
case of SN 2009ib we obtained tpT — to = 141±2 days, which 
is about 8 — 33 days longer than those of the comparison SNe. 

2.1 Reddening and Colour Curves 

The galactic extinction in the directi on of SN 2009ib is low, 
E(B — U)mw = 0.0257± 0.0002 mag (ISchlaflv fc Finkbeined 
I2OIII) . Using the classification spectrum taken with the 
Gemini South Telescope (see Sect. [Sj, we measured the 
equivalent width of the Na l D line at the redshift of the 
host galaxy as EW (NaD) = 0.838 ± 0.018 A (the compo¬ 
nents of the multiplet are not resolved in the spectrum). 
At the redshift of the Milky Way Na l D absorption can¬ 
not be detected in the spectrum, i n agreement with a low 
E{B — V)mw- Using the equation of IPoznanski et ^ ll2012l) . 
the equivalent width of Nal D leads to the host galaxy red¬ 
dening of E{B — U)host = 0.131 ± 0.025 mag. Together with 
the Milky Way component we estimate the total reddening 
as E{B — U)tot = 0.16 ± 0.03 mag. Note, however, that in 
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Figure 2. Optical and near-infrared light curve of SN 2009ib. The phase is given relative to the estimated explosion epoch, to = 
55041.3 MJD (Sect l4.1.1ll . 


several cases the Na l D feature was shown to be an unreli¬ 
able indicatior of e xtinction (see e.g. IPoznanski et al.ll201ll : 
iPhillips et al.llioi^ l. 

Fig. 0] shows the {B - V)o, (V - R)o, (V - I)o, {V - 
J)o, and {V — H)o colours of SN 2009ib, corrected with 
the reddening above, and compared with the colour curves 
of other SNe II-P. The colours of SN 2009ib do not differ 
significantly from those of other II-P SNe. 

During the second half of the plateau phase the colours 
of the SNe in our sample remain roughly constant, but when 
this phase ends they become signihcantly redder. The colour 
of SN 2009ib is mostly constant about 10 — 30 days longer 
than those of the other SNe and only then becomes as red¬ 
der, which is a manifestation of the longer plateau phase. 
During the nebular phase the colours are similar to those of 
SN 2004et. 


3 SPECTROSCOPY 

Optical spectroscopic observations of SN 2009ib were taken 
with the Gemini South Telescope (GST+GMOS), the New 
Technology Telescope (NTT-I-EFOSC2) and the Very Large 
Telescope (VLT-I-FORS2) between days 11 and 263 after 
explosion (for the summary of the observations see Table 
[2J. The image reduction was carried out using standard 
IRAF tasks. After the basic reduction steps the spectrum 
was extracted and wavelength calibrated using comparison 
lamp spectra. We used spectra of standard stars observed 


Table 1. Explosion epoch, reddening, distance and plateau 
length of the SNe used for comparison throughout the paper. 


SN 

^0 

(MJD) 

E{B - V) 
(mag) 

D 

(Mpc) 

^0 ^PT 

(days) 

Ref. 

1999em 

51476.5 (1.0) 

0.10 (0.05) 

11.7 (1.0) 

112 (3) 

1,2,3 

1999gi 

51517.8 (3.0) 

0.21 (0.10) 

11.1 (2.0) 

128 (5) 

4 

2002gd 

52551.5 (2.0) 

0.06 (0.01) 

37.5 (3.5) 

113 (7) 

5,6 

2003Z 

52664.5 (4.5) 

0.04 (0.01) 

21.9 J.6) 

133 (8) 

5 

2004et 

53270.0 (0.2) 

0.41 (0.01) 

4.8 (0.4) 

123 d) 

7,8 

2005CS 

53548.5 (0.5) 

0.05 (0.01) 

8.4 (0.7) 

125 (6) 

9,10 

2008in 

54825.1 (0.8) 

0.10 (0.10) 

12.2 (1.9) 

108 (3) 

11 

2009N 

54847.6 J.2) 

0.13 (0.02) 

21.6 J.l) 

109 (2) 

12 

2009js 

55115.4 (5.5) 

0.36 (0.12) 

21.7 (1.8) 

117 (7) 

13 


Referen ces: fll |Hamuy et al.| <l200l|h (2) [Leonard et al.l 

ll2002a^. Odbeo nard et alj ||200^. (4) [ L eonard et alj ll2002bl') . 
(5) ^iro^^ (6) I Poznansk^^al] |200^, 

Takdts J2^^, (8) iMa guir^^^^T <l2010li . (9) 

Pastorelloe^^ 11200^. (10) Vink^^t al.l <[201^, (H) 
^P ^^^alJ 1 I 2 OII ), fl2) lT^akats et al.lll2014ft . T3) [Gandhi et al.l 


J2013h . 


tp'T is the epoch of the middle of the transition from the 
plateau to the tail. 


on the same night in order to flux calibrate the SN spec¬ 
trum. The flux was then checked against the photometric 
measurements taken at a close epoch and corrected using a 
scaling factor if necessary. The spectra are shown in Fig. [5] 
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Figure 4. The reddening corrected colour curves of SN 2009ib in 
comparison with those of other SNe II-P. 


and[6l In the case of the first spectrum taken with GST there 
was no flux standard star observed, therefore we used the 
standard taken during the night of our second GST obser¬ 
vation, on August 14, for flux calibration. The continuum 
of both spectra, indicate unreasonably high temperatures, 
in disagreement with the photometric measurements. Since 
we have another spectrum taken by NTT at the same epoch 
as the second GST spectrum, we adjusted the continuum of 
both of the GST spectra with the help of this NTT spec¬ 
trum, resulting in more realistic but quite uncertain continua 
(Fig. El). 

The classification spectr um of SN 2009ib wa s taken on 
Aug. 9.1 UT with the GST llPignata et al.|[2009l l. Lines of 
the Balmer-series of Hi and He l A5876 are present. Ac- 
cording to the SNIP program (Supernova Identification, 
iBlondin fc Tonr^l2007h the spectrum is similar to that of 
SN 2005cs about a week after explosion. By day 16 metal¬ 
lic lines had appeared such as Fell, Tin, Sin, Can, while 
He I was not visible any more. The third spectrum was taken 
at the phase of 46 days, when the metallic lines were much 
stronger and lines of Nal D, Sc ll, Ball and O l also appeared. 
During the rest of the plateau phase these lines persisted, 
and also remained present after the luminosity drop, at the 


> 



Figure 7. Comparison of the expansion velocities of SN 2009ib 
measured from the Fe II A5169 to those of the normal SNe II-P 
1999em, 1999gi and 2004et, the subluminous SNe 2002gd, 2005cs, 
2003Z and the intermediate luminosity SNe 2008in and 2009N. 
The insert shows Doppler-velocities from the absorption minima 
of H«, H/3, Fell A5169 and Sen A6245 lines of SN 2009ib. 


phase of 263 days. In our last two spectra, taken on days 
219 and 263, nebular lines are also detected. The emission 
component of the Call IR triplet became prominent, and for¬ 
bidden lines, such as [Call] AA7291,7323, [Oi] AA6300,6364, 
[Fell] AA7155,7172 are present (FigEJ- 

The spectra of SN 2009ib have narrow lines with low 
velocities. We measured the Doppler-velocities of Hq, H/3, 
Fell A5169 and Sen A6245 lines during the photospheric 
phase by fitting a gaussian to their absorption minima. The 
velocity curves are shown in the insert of Fig. [T] We compare 
the velocities of SN 2009ib to those of other SNe in Fig[3 
measured from the absorption minimum of the Fell A5169 
line. The velocity curves of these objects have a high scatter, 
the “normal” SNe II-P, like SNe 1999em, 1999gi and 2004et 
have higher velocities than the rest of the sample. The ve¬ 
locities of SN 2009ib are somewhat higher than those of the 
intermediate luminosity SNe 2008in and 2009N. 

In Fig. E] we compare the -1-83 day spectrum of 
SN 2009ib with the spectra of normal and subluminous 
SNe taken at similar phases. The shape, width and rela¬ 
tive strength of the features in the spectrum of SN 2009ib, 
especially the lines of Hi, Ball and Fell, appear to be the 
most similar to those of the subluminous SN 2002gd and 
differ somewhat from those of SN 2009N, and from those of 
the other subluminous SN included in this comparison. It 
demonstrates that the observational properties of SNe II-P 
are diverse, they do not form clearly distinguishable groups. 


3.1 NIR spectra 

Two NIR spectra were taken with NTT (-I-SOFI), on days 
-1-48 and -1-84 (Table [2)). The reduction steps were carried 
out using standard iraf packages. Each night several pairs 
of spectra were taken at different positions along the slit. 
These pairs were subtracted from each other in order to 
remove the sky background, the subtracted images were then 
added together. The SN spectrum was extracted from the 
co-added image. The wavelength calibration was carried out 
using arc lamp spectra that was taken before and after the 
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Table 2. Summary of the spectroscopic observations. 


Date 

MJD 

Phase^ 

(days) 

Instrument set-up 

Wavelength range 

(A) 

09/08/2009 

55052.4 

11.1 

GST + GMOS + R400 

3900-8100 

13/08/2009 

55057.3 

16.0 

NTT + EFOSC2 + gr#ll + gr#16 

3500-10100 

13/08/2014 

55057.4 

16.1 

GST + GMOS + R400 

3900-8100 

12/09/2009 

55087.4 

46.1 

NTT + EFOSC2 + gr#ll + gr#16 

3500-10100 

14/09/2009 

55089.3 

48.0 

NTT + SOFI + GB + GR 

9500-25200 

19/10/2009 

55124.3 

83.0 

NTT + EFOSC2 + gr#ll+gr#16 

3500-10100 

20/10/2009 

55125.3 

84.0 

NTT + SOFI + GB 

9500-16400 

22/11/2009 

55157.2 

115.9 

NTT + EFOSC2 + gr#ll + gr#16 

3500-10100 

24/12/2009 

55190.2 

148.9 

NTT + EFOSC2 + gr#ll + gr#16 

3500-10100 

24/01/2010 

55221.2 

179.9 

NTT + EFOSG2 + gr#ll 

3500-9500 

04/03/2010 

55260.1 

218.8 

VLT + FORS2 + 300V + 3001 

3500-11000 

17/04/2010 

55304.0 

262.7 

NTT + EFOSG2 + gr#13 

3500-9500 


“ relative to the estimated date of explosion, to = 55041.3 MJD 



Rest wavelength (Angstrom) 


Figure 5. The spectral evolution of SN 2009ib during its photosheric phase. Note, that in the case of the Gemini spectra the flux 
calibration is approximate (see text). In the figure we identify the strongest features present in the spectra. The wavelength positions of 
telluric features are marked with © symbol. 


SN spectra. The strong telluric features were removed with 
the help of the spectrum of a G-type telluric standard star 
observed close in time, with the same instrumental setup and 
at similar airmass as the SN. We also used the spectrum 
of the telluric standard for the flux calibration of the SN 
spectra that were later checked against the NIR photometry 
from the nearest epoch and corrected when necessary. 

The NIR spectra of SN 2009ib are shown in Fig. We 
marked the position of the spectral lines commonly seen in 
NIR spectra of SNe II-P. Along with the H i features, lines of 


Sr II, Cl and Fell are visible. The first NIR spectrum covers 
also the range redwards from 16000 A, where we marked the 
position of Hel A20580, however the line is not visible. In 
Fig. Ho] we compare the NIR spectrum of SN 2009ib taken at 
phase +84 d ays to those of SN 2009N and the subluminous 
SN 2009md llFraser et al.ll201ll ) at similar phases. While the 
Sr II lines seem to be weaker, the Cl A10691 feature in the 
spectrum of SN 2009ib is significantly stronger and broader 
than those of the other two. The shapes of the emission 
components of the Hi lines of SNe 2009N and 2009ib are 
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Figure 6. The spectra of SN 2009ib during its tail phase. The lines that are typically present in the nebular phase spectra of SNe II-P 
are marked. A spectrum of the subluminous SN 2005cs at similar phase is also shown for comparison. 


similar, but somewhat wider than those of SN 2009md. The 
feature at 10550 A (whi ch was present in the spectrum of 
SNe 2009N and 2008in, IXakats et al.l[20n i also appears in 
the spectra of SN 2009ib (marked with an arrow in Fig llOII . 
however less pronounced, maybe due to the st rong, wide C l 
featur e next to it. This line was identified bv lTakats et al.l 
l|2014h as either a high-velocity component of Hel A10830 or 
a Si I feature, however they were unable to draw a definite 
conclusion (see their Fig.16). 


4 DISTANCE 

The distance to the host galaxy of SN 2009ib, NGC 1559, has 
been measured several times. The results, however, show a 
high scatter, falling in the range between 12.6 and 22.0 Mpc 
(see Sect. I4.4|l . 

There are two methods that allow us to measure dis - 
tances of S Ne II-P directly, EPM dKirshner &: Kwanlll974l ~l 
and SCM jHamuv fc Pintal2002ll . We apply both methods 
on the data of SN 2009ib (Sect l4.l1 and l4.2Tl and EPM on the 
data of SN 1986L, an earlier SN that exploded in the same 
galaxy (Sect 14.3|) . 

4.1 EPM 

EPM is a variant of the Baade-Wesselink method, its appli¬ 
cation requires photometric and spectroscopic monitoring 
throughout the first half of the plateau phase. A great ad- 
vantange of this method is that it does not require calibra- 



Rest Wavelength (Angstrom) 


Figure 10. Comparison of the -|-48d NIR spectrum of SN 2009ib 
(red) to those of SNe 2009N and 2009md at similar phases. The 
po sition of an in t erestin g feature at 10550 A, that was identified 
by iTakats et al.l ||2014| ') as either a high-velocity component of 
He I A10830 or a Si l feature is marked with an arrow. 

tion using objects with known distances. Assuming that the 
optically thick ejecta expands homologously and radiates as 
a diluted blackbody, the method relates the apparent radius 
of the photosphere to its physical radius. The method has 
been discussed extensively in the literature, so for details we 
refer the reader to the relevant papers (e.g. iHamu 3[20i|; 
[Leonard et al.|[2002al : iDessart fc Hillie^l20()^ ). 

To calculate the physical radius of the photosphere we 
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Figure 8. Comparison of the spectra of selected SNe II-P at similar phases. The spectra are sorted by the brightness of the SNe, the 
brightest one being on the top, and includes those of the normal SNe 1999em, 2004et, the subluminous SNe 2002gd, 2003Z, 2005cs and 
the intermediate luminosity SN 2009N. 



Figure 9. The near-infrared spectra of SN 2009ib, taken 48 and 84 days after explosion. 
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Table 3. The calculated parameters of SN 2009ib needed to apply 
EPM: the apparent radius (0), photospheric temperature (T), di¬ 
lution factor (C(^))i the photospheric velocity (u), measured 
from the absorption minimum of Fe II A5169 and interpolated to 
the epochs of the photometry. 


MJD 

t-to“ 

1-1-2 

(days) 

e 

(10® i^) 

T 

(K) 

C(F) 

V 

(km/s) 

55054.3 

13.0 

3.481 

10820 

0.52 

6999 

55057.3 

15.9 

3.959 

9278 

0.54 

6181 

55063.3 

21.9 

4.860 

7578 

0.62 

5107 

55064.3 

23.0 

4.656 

7622 

0.61 

4965 

55066.3 

24.9 

5.172 

6780 

0.68 

4729 

55066.3 

25.0 

4.895 

7097 

0.65 

4729 

55067.3 

26.0 

4.907 

7110 

0.65 

4621 

55068.3 

26.9 

4.941 

6956 

0.66 

4522 

55069.3 

27.9 

5.538 

6199 

0.75 

4433 

55070.3 

28.9 

4.934 

6867 

0.67 

4339 

55072.3 

30.9 

5.091 

6497 

0.71 

4176 

55072.4 

31.1 

5.121 

6579 

0.70 

4172 

55074.3 

32.9 

5.239 

6205 

0.75 

4028 

55075.4 

34.0 

5.124 

6540 

0.71 

3956 

55077.3 

35.9 

5.398 

6022 

0.77 

3832 

55079.3 

37.9 

5.864 

5512 

0.87 

3714 


“ relative to the estimated date of explosion, to = 
55041.3 MJD 


used equation (1) of lHamuv et al.l |2003), and applied it to 
the BVI band photometry. We adopt the velocities mea¬ 
sured from the absorption minimum of the Fell A5169 fea¬ 
ture as the photospheric velocities. We interpolated them to 
the epochs of t he photometry us i ng th e relation given in the 
equation (2) of iTakats fc Vin^ (l2012l l. The dilution factors 
can be determined as a function of the photospheric tem¬ 
perature from atmospheric model cal culations. So fa r there 
hav e been two such studies, those of lEastman et ahl (Il996l l 
and iDessart fc Hillierl ll2005h. In this paper we adopt ed the 
more recent value s, those of Dessart &: Hillien (l2005ll . Note 
that the values of lEastman et ahT i 1996ll are systematically 
lower, leading to shorter distances. The quantities calculated 
to apply EPM can be found in Table [S] The resulting dis¬ 
tance is Depm = 21.3 ± 2.2 Mpc, while the explosion date 
determined via EPM is to = 55041.3 ± 3.1 MJD. 


4-1.1 Confirming the date of explosion 

In order to test the explosion epoch obtained via EPM, we 
estimated this date by matching the plateau phase spectra 
of SN 2009ib to those of other si milar objects. We used th e 
applications snid and GELATcjfl (iHarutvunvan et al.ll200^ . 
both of which compare the given spectrum to extended 
databases of SN spectra. 

First, we selected a sample of comparison SNe that 
have good spectral covarage and their explosion epochs are 
well-constrained. Then we narrowed this sample by com¬ 
paring their spectra to those of SN 2009ib, in order to fil¬ 
ter out those objects that have significantly different spec¬ 
tral evolution. The final sample of templates contains SNe 
1999em, 2002gd, 2004et, and 2009N for snid, and SNe 


® https;//gelato.tng.lac.es 


1999em, 2004et, 2012A and 2007od in the case of GELATO. 
We matched each plateau phase spectra of SN 2009ib to the 
templates of these comparison SNe one by one, and calcu¬ 
lated the explosion epoch taking into account the quality 
of the fits, i.e. the rlap parameter in the case of snid and 
the Quality of Fit parameter in the case of GELATO. The 
systematic errors are calculated as the RMS of the explo¬ 
sion epochs from the different template SNe. In the case 
of SNID the mean value of the determined explosion epochs 
was tsnid = 55044.7 ± 6.5 (random) ± 1.3 (systematic). With 
GELATO the systematic error of the results is higher, we ob¬ 
tained tgeiato = 55040.3±7.7 (random) ±3.2 (systematic). In 
spite of the high uncertainties due to the slow spectral evo¬ 
lution of SNe II-P during the plateau phase, these results 
confirm the epoch obtained via EPM, which therefore we 
adopt as the explosion epoch of SN 2009ib. 


4.2 SCM 


Another method to calculate distances to SNe II-P is SCM 
which is based on an empirical correlation between the 
brightness of the SN and it s expansion velocity in the 
middle of the plateau phase dHamuv &: Pintd l2002tl . The 
method requires calibration using objects with well-known 
distanc es. The original version has been revised multiple 


times iHamuyl [2005 


Nugent et al.l l2006l : IPoznanski et al.l 


I 2 OO 9 I : lOlivares et al.l Her^ we apply s everal versions 


simi larly as was done in Takats et al.l (l2014l l. 

IPoznanski et al.l 1 I 2 OO 9 II calibrated the equation 

Mi- a- log ) ± Ri{{V - I) - {V - I)o) - mi 

= -5 • \og{HoD) 


and defined Mi = —1.615 ± 0.08 mag, a = 4.4 ± 0.6, Ri = 
0.8 ±0.3 and (V — I)o = 0.53 mag. In the case of SN 2009ib, 
we measured m/ = 15.368 ± 0.020 mag, {V — I) = 0.800 ± 
0.032 mag, and estimated iiFe(50d) = 3247 ± 200 km/s. As¬ 
suming a Hubble constant Ho = 72 km s“^ Mpc”*^, we cal- 
cula ted the distance as D scm.i = 21.48 ± 2.00 Mpc. 

lOlivares et al.l 1 I 2 OIOII examined the dat a of 37 nearby II- 
P SNe , and applied the same expression as iHamuv fc Piiit^ 
II 2 OO 2 D : 

m ± a log(uFe/5000) — l3{V — I) = 5 log HoD ± zp, (2) 


but using the magnitudes and velocities measured 30 days 
before the middle of the transition phase (Ipt). They cal¬ 
ibrated the formula using differe nt bands, the values of fi 
and zp can be found in Table 6 of lOlivares et al.l (I 2 OI 0 I I. In 
Section [ 2 ] we determined tpT = 55182.7 MJD (141 days af¬ 
ter explosion). At 30 days before this date we measured the 
values v(llld) = 2010 ±300 km/s, ms = 18.49 ± 0.06 mag, 
mv = 17.10 ± 0.03 mag, mu = 16.39 ± 0.02 mag and 
mi = 15.99 ± 0.02 mag. These values led to the distances 
23.70 ±3.50 Mpc, 26.74 ±3.29 Mpc and 28.09 ±2.93 Mpc in 
B, V and I bands, respectively. The mean of these results 
is Dscm,2 = 26.19 ± 3.25 Mpc. 

We also apply the version of iMaguire et al. I ll2010ll . 
who extended the te chnique to the N I R ran ge. They used 
the same formula as IPoznanski et (l2009l l (equation [T]). 
Their calibration in J band yielded the val ues ol M. 7 = 
—2.532 ± 0.250 mag and a = 6.33 ± 1.20. [Maguire et al.l 
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ll2010fl ch ose to calculate with th e value of Rv = 1-5, ob¬ 
tained by IPoznanski et alJ (120091 '). instead of using it as a 
fitting parameter. In the case of SN 2009ib, we measured 
mj = 14.88 ± 0.16 mag, {V — J) = 1.29 ± 0.16 mag, 
and used UFe(50d) = 3247 ± 200 km/s at 50 days af¬ 
ter the explosion. In this way we obtained the distance 
15scm,3 = 23.16 rh 3.48 Mpc. 

The three variants of SCM provided slightly different 
distances. We adopt their mean, Dscm = 23.6 ± 1.7 Mpc as 
the SCM distance to SN 2009ib. 


4.3 Distance to SN 1986L 


SN 1986L was another type II-P supernova that exploded 
in NGC 1559. It was discovered at a very early phase and 
its dat a were also exploited to estimate the host galaxy dis¬ 
tance. ISchmidt et al.l l| 19941 ') using an early implementation 
of the EPM and the dilution factors of lEastman et all (|l996l 'l 
(which were at that tim e not ye t publi shed) estimated the 
distance as 1 6 ± 2 Mpc. Harnuyl |2001^ also used the dilu¬ 
tion factors of lEastman et all ( 1996l l and a cross-correlation 
technique for measuring v elocities and o btained the distance 
of 11.2 Mpc. A c cordi ng to lHamir^ 1 I 2 OOII I. the higher value of 
ISchmidt et al.l d 19941 ') is due to the shallower velocity curve 
they adopted. 

We used photometric measurements in BV bands (M. 
Phill ips, private communication; see also lAnderson et al.l 
l2014h . and velocities measured from the Doppler-sh ift of 
the Fell A5169 line obtained from ISchmidt et al.l (Il994l l and 
iHamuvl (1200 ih . We interpolated the velocit ies to the epochs 
of the photome t ry us ing the relations of iTakats &: Vinkg 
(I 2 OI 2 II . Tllamuvl (1 200 il l found that the host galaxy extinc¬ 
tion is negligible, therefore we only take into account the 
galactic extinction in the direction of SN 1986L, which is 
E(B — V)mw = 0.0258±0.0003 mag l|Schlaflv fc Finkbeinerl 

M). 


We applied EPM the same way a s in the case of 
SN 20 09ib, using the dilution factors of IPessart fc Hillierl 
l|2005h and photometry in BV band. The calculated param¬ 
eters can be found in Table ID We obtained the distance as 
Dsen = 17.2 ± 3.0 M pc, l onger than the previous results of 
ISchmidt et al.l (|l994l l and IHamuvl l|200 j l. but shorter than 
the distance determined from the data of SN 2009ib. The 
difference from the earlier results is mainly due to the di¬ 
lution factors, which are based on atmospheric models that 
greatly improved since the first applications of EPM. The 
lack of I band data usually leads to more uncertain tem- 
perature estimate s and slightly different distances (see e.g. 

I Jones et al.ll200^ 1. In the case of SN 2009ib, taking into ac¬ 
count only the BV bands would lead to a distance only 4 
per cent lower than with BVI bands. 


4.4 Average Distance 

As mentioned in the beginning of this Section, there have 
been several previous distance measurements to the host 
galaxy of SN 2009ib, and those results have a high scatter. 

Mos t of the measuremen ts used the Tully-Fisher (TF) 
method (iTullv fc Fisheriri977l l. Since the TF method was ap¬ 
plied by only a few groups, several times either by recalibrat¬ 
ing the method or adding newer data, we decided to select 


Table 4. The calculated parameters of SN 1986L needed to apply 
EPM: the apparent radius (6), photospheric temperature (T), di¬ 
lution factor (C(T))i and the photospheric velocity (v), measured 
from the a bsorptio n min imum of Fell A5169 by I Schmidt et ^ 
(Il994h and|H amuvl (l200lh and interpolated to the epochs of the 
photometry. 


MJD 

t-to ^ 

1 + 2 
(days) 

e 

fins Ahl-) 

Mpc-' 

T 

(K) 

C(T) 

V 

(km/s) 

46712.4 

1.8 

3.495 

25920 

0.42 

30154 

46714.3 

3.7 

6.569 

13715 

0.46 

18863 

46715.3 

4.7 

6.537 

14171 

0.46 

16325 

46716.2 

5.6 

7.549 

11944 

0.49 

14720 

46717.3 

6.7 

7.495 

12084 

0.49 

13253 

46729.3 

18.6 

8.516 

8615 

0.62 

7404 

46730.4 

19.7 

8.993 

7770 

0.69 

7173 

46732.4 

21.7 

9.289 

6839 

0.80 

6798 

46735.2 

24.5 

8.713 

7097 

0.76 

6357 

46736.2 

25.5 

9.954 

5894 

0.98 

6218 

46737.4 

26.7 

9.484 

6651 

0.83 

6062 

46738.4 

27.7 

9.428 

6150 

0.92 

5940 

46740.3 

29.6 

9.223 

5988 

0.96 

5726 

46742.3 

31.6 

9.072 

5711 

1.03 

5523 

46744.3 

33.6 

10.447 

5082 

1.24 

5340 

46748.2 

37.4 

9.836 

5090 

1.23 

5027 

46750.3 

39.5 

10.535 

4653 

1.43 

4878 

46754.3 

43.4 

10.832 

4517 

1.51 

4633 

46756.3 

45.4 

10.699 

4440 

1.55 

4519 

^ relative to the epoch of discovery, to = 

46710.6 MJD 

llEvans et alJll986h. 




only the latest paper from each group (IWillick et al.lll997l: 


ISpringob et al.ll2007l : ITuUv fc Courtoisll2012l '). and then av¬ 
erage their results. In this way we estimated the TF distance 
as 15.7 ± 1.5 Mpc. 

The distance to this galaxy was also determined using 
previous SNe. As discussed in Sect. 14.31 the EPM method 
was applied to the Typ e II-P SN 1986L bv ISchmidt et ahl 
1 I 1994 II and by [h amuvl ll200lf) . but after reanalysing these 
data we obtained a distance somewhat higher than the pub¬ 
lished values. The MLCS2k2 method was used for the Type 
la SN 2005df leading to a distance of 22.0±2.0 Mpc (?). The 
redshift distance is 14.9±1.0 Mpc, however it has significant 
uncertainties due to the effect of peculiar motions, therefore 
we exclude it from our analysis. Considering the TF dis¬ 
tance, that of SN 2005df, our EPM result for SN 1986L, 
and the distance to SN 2009ib determined with EPM and 
SCM (Table [5} , we conclude the mean of these five results, 
D — 19.8 ± 2.8 Mpc as the distance to NCC 1559. 


5 BOLOMETRIC LUMINOSITY AND NI MASS 

Using the distance determined in Sect. |4l we can calculate 
the bolometric luminosity of SN 2009ib. Between days 13 
and 104 we have observations covering the wavelength range 
from u' band up to H band (see Sect. [2}. Using low order 
polynomials, we interpolated the magnitudes in u'BRIJH 
bands to the epochs of the observations in V band. After 
correcting for the reddening, we converted these magnitudes 
to fluxes, and integrated them using Simpson’s rule. The in¬ 
frared range redwards from H band was estimated assuming 
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Table 5. The distances to the host galaxy of SN 2009ib obtained 
from the literature as well as those determined in this paper. The 
mean of the values is also shown. 


Method 

D 

h 

References^ 


(Mpc) 

(mag) 



Tully-Fisher 

15.7 

(1.5) 

30.72 

(0.21) 

1,2,3 

SN la (2005df) 

22.0 

(2.0) 

31.72 

(0.20) 

4 

EPM (1986L) 

17.2 

(3.0) 

31.18 

(0.37) 

this paper 

EPM (2009ib) 

21.3 

(2.2) 

31.64 

(0.21) 

this paper 

SCM (2009ib) 

23.6 

(1.7) 

31.77 

(0.15) 

this paper 

Mean 

19.8 

(3.0) 

31.48 

(0.31) 



^ (1) IWillick et al.l ifl^h ^2^1 ISoringob et al.l J20Q7ri . (3) 
ITuIIv ^^^our^isl ||201 j) . (4) iMilne et al.l l|2010lV 

that it follows the Rayleigh-Jeans tail of a black body spec¬ 
trum. After day 104 we lack observat ions in the NIR bands . 
At this phase we used the formulae of lMaeuire et al.l ll2010ll , 
who estimated the bolometric luminosity from the R and 
V magnitudes. We find that during the plateau phase this 
estimation agrees with the luminosity that we calculated, 
therefore we use them at the later phase, after day 104. 
The obtained uvoir bolometric luminosity curve is shown in 
Fig. mi 

Since the observational coverage for some of the SNe is 
incomplete, for comparison purposes we computed the lumi¬ 
nosities by integrating only the BVRI fluxes. The distances 
of the comparison SNe can be found in Table □ Fig. mi 
shows the comparison of the quasi-bolometric luminosity of 
SN 2009ib to those of other SNe II-P. In this comparison 
we included several SNe that were categorized as sublumi- 
nous (SNe 2002gd, 2005cs) as well as “normal” events (SNe 
1999em, 1999gi, 2004et) and “intermediate” ones, such as 
SNe 2008in, 2009N, 2009js. The luminosity of SN 2009ib 
during its plateau phase is similar to that of SN 2009N. 
Authors, w ho examined a larger number of light cur ves of 
SNe II (e.g. I Anderson et al.ll2014l : lenders et al.ll201^ have 
found that the luminosity distribution of SNe II-P is contin¬ 
uous, which is reflected in our Figure. 

Fig. Illl also shows that t he plateau phase of SN 2009ib 
is the longest in the sample. lAnderson et al.l (120141 ') found 
a weak anticorrelation between the plateau length and the 
absolute brightness including only 15 SNe in the sample. It 
is clear, that SN 2009ib does not fit in this trend. 

The ®®Ni mass produced in the explosion can be esti- 
mated f r om th e tail phase luminosity. We used the relation of 
iHamuvl (120031 ') and obtained the value of 0.046 ± 0.015 Mq. 
This value is signihcantly higher than those of the sublumi- 
nous and intermediate SNe II-P ISec. 16.41) . 


6 PROGENITOR PARAMETERS 

Even though a large number of SNe II have been discovered 
and studied, we still do not fully understand the exact na¬ 
ture of their progenitors and explosions. In a number of cases 
direct identification of the progenitor star on archival pre¬ 
explosion images have revealed that SNe II-P are originated 
from red sup ergiant stars in the mass range 8 — 16 Mq (eg. 
ISmarttll2OT9 1. In other cases, the parameters of the progen- 
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Figure 11. Comparison of the evolution of the BVRI quasi- 
bolometric light curve of SN 2009ib (black filled circles) to those 
of other SNe II-P. We also show the uvoir bolometric curve of 
SN 2009ib (black open circles). 


itor and its explosion have be en estimated via hydrodynam- 


ical model calculations (e.g. Litvinova&^Jadj;ozhW 


Zamgieri^j^^ ]j2QQ^; IPumo fc ZamDierill201ll : IUtrobiril l 


Dessart et al.l 201fll ). 


1983; 


20071 : 


In this section we investigate the properties of the pro¬ 
genitor of SN 2009ib first by examining the archival pre¬ 
explosion images of its location taken by the HST, then, 
independently, by hydrodynamical modelling of the SN ob¬ 
servables and, finally, we compare the light curve and spec¬ 
tra of SN 2009ib to th e non-LTE time-depe ndent radiative- 
transfer simulations of lPessart et al.l (l2013l ). 


6.1 HST images 

We downloaded the pre-explosion images of the SN site from 
the HST archival. These images were taken with WFPC2 
using the filters F45QW, F&Q&W and ESldlF with the to¬ 
tal exposure time of 320 s each, and the pixel scale of 
0.0996” pixel”^. The images were obtained on 2001 August 
2 (Proposal ID 9042, PI. S. Smart!). 

In order to measure the precise position of the super¬ 
nova, we used post-explosion images taken with VLT-|-NaCo 
in Ks band (PI: Smart!), on 13 August 2009, about 16 days 
after explosion, downloaded from the ESO archiv^. These 
images have the pixel scale of 0.05”pixel“^, the field of view 
of 56” X 56”, and the total exposure time of 3960 s. We 
carried out dark- and flat-held corrections on all 22 images 
using the standard tasks of iraf. The sky was substracted 
by using off-source sky images taken before and after the on- 


^ https://archive.stsci.edu 

® http://archive.eso.org/eso/eso_archive_main.html 
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source images. The sky-substracted frames were then aligned 
and combined. 

In order to achive precise relative astrometry, we geo¬ 
metrically transformed the HST images to the NaCo image 
using the GEOMAP and GEOTRAN tasks of iraf, measuring 
the positions of 11 point sources that are present on both 
images. We measure the SN position in the NaCo image 
with the IRAF tasks daofind and imexamine. The error of 
the position of the SN was calculated from the uncertain¬ 
ties of both the geometric transformations and the position 
measurements as 160 mas. 

We used the dolphot prograrr0 to detect and measure 
sources present in the pre-explosion images. Dolphot de¬ 
tected two sources close to the SN position. Both sources 
were detected at > 5cr with the filters F606W and FSldlT, 
but not with filter F4501T. 

Fig |12l shows the post-explosion image, as well as the 
same field in the pre-explosion image (in FSIAW band). 
The third image on this figure zooms in to the surround¬ 
ings of the SN position, with the detected sources marked 
as “A” and “B”. This figure shows, that source “A” is a pos¬ 
sible progenitor of the SN. Table [6] contains the positions of 
these sources relative to the SN, as well as their brightness 
measured by dolphot. We converted these magnitudes to 
the Johnson-C ousins system using the constants provided by 
iDolphirJ (l2009ll and corrected for the reddening fSect. ETT1) . 
The obtained VI values can also be found in Table |6l 

The spectral type of these sources (assuming they 
are single stars) wer e esimated from t h eir V — I colours 
using Table 4 o f Maund &: SmarttI (l2005ll (see also 
iDrilling fc Landoltll20od l. We found that both sources are 
likely type G stars. Therefore if source A is the progeni¬ 
tor, it is more yellow, than the the RSG stars that are the 
usual prog enitors of SNe II-P . We u sed the bolometric cor¬ 
rections of iMaund fc Smarttj (l2005ll to calculate the lumi¬ 
nosity of source A. We obtained the value of \og[L/L q) = 
5.04 ± 0.2 dex. 

lAnderson et al.l il2010l l examined the site of SN 1986L, 
which exploded in the same galaxy as SN 2009ib, and mea¬ 
sured the metallicity as (12 -|- log(0/H)) = 8.51 ± 0.13, a 
solar-like metallicity. Using the spectroscopy image taken 
with VLT (Sect[31), we extracted the spectrum of a Hll 
region close to SN 2009ib. The distance of this region to 
SN 2009ib is 17.1”, which corresponds to 1.6 kpc. It is 
closer than SN 1986L, the distance of which is 65.5”. From 
the spectrum of the Hll region we measured the ratios of 
[Nii]A6583/Ha (N2) and ( fO iii]A5007/Hm/N2 (0 3N2), and 
following the pr ocedure oflAnd erson et al.l (l2010l l and using 
the relations of IPettini &: Pagell ( 2004 1 we obtain the N2 
metallicity (12 -|- log(0/H))N2 = 8.44 ± 0.21 and the 03N2 
metallicity (12 + lo g(0/H))o3N2 = 8.73 ± 0.17. This value 
is similar to that of lAnderson et ^ (l2010h and close to the 
solar metallicity. 

In Fig. [13] we plotted the position of source A on a 
Herzsprung-Russel (HR) diagram toge ther with stellar evo- 
lution tracks from the STARS code fsee lEldridge et ^l2008l . 
and the references therein) with solar metallicity. This dia- 


® DOLPHOT is a stellar photo metry pa c kage t hat was adapted from 
HSTphot for general use by IDoIdIm ||200CI| ~). 


Table 6 . Position and brightness of the sources found by 
DOLPHOT in the pre-explosion HST images close to the SN po¬ 
sition. 



Source A 

Source B 

Aa/A(S (mas)“ 

4/50 

20/310 

F606W (mag) 

23.91 (0.11) 

24.05 (0.17) 

F814W (mag) 

23.20 (0.12) 

23.12 (0.13) 

V^ (mag) 

23.79 (0.18) 

24.02 (0.24) 

F (mag) 

22.94 (0.17) 

22.83 (0.21) 

Spectral type‘s 

GO-Gl (F6-G6) 

G8-K0 (G1-K3) 


“ distance to the SN position. 

^ reddening corrected values. 

in parenthesis we give the spectral type range estimated 
considering the uncertainties of the VI magnitudes. 


gram shows that source A is on the evolution track of the 
20 Mq star. 

In a few cases there have been a yello w source found 
at the SN postition in pre-explosion images. iLi et al.l (I 2005 II 
found the progenitor of SN 2004et t o be a yellow superg iant 
with a ZAMS mass of 15^2 ^^01 but lSma rtt et all (|2009|) ar¬ 
gued that that source is not a single star. Elias-Rosa et al.l 
1 I 2 OO 9 II also found a yellow source at the position of the pro¬ 
genitor of SN 2008cn. Other than having a yellow supergiant 
with a ZAMS mass of 15 ± 2 Mq as a progenitor, they offer 
alternative explanations for the yellow colour of the source. 
They consider the possibility that it is not a single star, 
instead a blending of two or more star. Examining model 
stellar spectra, they showed that the yellow colour can re¬ 
sult e.g. from a blend of a type BO with Mini ~ 40 Mq and 
a type Ml star with Mini ~ 15 Mq. In fact, by taking very 
late time images of this SN with HST, iMaund et al.l (120151) 
showed that the pre-explosion yellow source contained a blue 
star that remained present after the explosion, and con¬ 
cluded that the progenitor was mos t likely an RSG star with 
Mini ~ 16 M^t^ . I n a s imilar case, iFraser et al.l (|201fll l and 
lElias-Rosa et ah lllw^) both arrived to the conclusion that 
the progenitor of SN 2009kr was a yellow supergiant star, 
with initial mass of 18 — 24 Mq and 15^^ Mm , respe ctively. 
Revisiting this SN at late times, iMaund et ah! (I 2 OI 5 II found 
that the yellow source is probably a small compact cluster, 
not an single source. 

In the case of SN 2009ib, we can examine a scenario sim¬ 
ilar to that of SN 2008cn. We can assume, that in the FUAW 
filter we detect the RSG progenitor, but due to its red color 
we do not see it in the FQQQW band, where the detection is 
due to a blue star, unrelated to the SN. Assuming that the 
RSG star has a temperature of ~ 3400 K, we determined the 
luminosity of the RSG star as log(L/LQ) = 5.12 ±0.14 dex, 
following the same procedure as before. We show the po¬ 
sition of this source on the HR diagram in Figure 1131 It 
shows that this RSG star would have had an initial mass of 

16±2Mq. 

We should also consider, that source A is not the pro¬ 
genitor, which, instead, is too faint to be visible in the pre¬ 
explosion images. Therefore we estimated the detection lim¬ 
its by running artificial star experiments using dolphot. 
We found that the upper limits at the position of the pro¬ 
genitor are 23.87 mag for EdhOlU, 24.30 mag for F606W 
and 23.25 mag for FSIAW. Correcting for the reddening 
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Figure 12. (a) The post-explosion image taken with VLT-hNaCo 15 days after the explosion, (b) A pre-explosion image taken by 
HST/WFPC2 with the filter F814W. The scale and orientation of the images (a) and (b) are the same. The circle on (b) shows the 
position of the SN with Icr uncertainty. The area marked with a square is enlarged in (c). The circle in (c) also shows the position of the 
progenitor, while A and B marks the sources detected with dolphot (see text and Table EJ. 


(Sect l2.ll and taking into account the distance to the SN 
(Sect[4|, these values correspond to the absolute magni¬ 
tudes Mf4,50w > —8.12 mag, Mpmew > —7.55 mag and 
Mf814w > —8.62 mag. 

We used the colo ur and bolometric corrections of 
iMaund fc SmarttI (l2005l l to convert the measured upper lim¬ 
its to bolometric magnitudes. Since in this scenario we have 
no precise colour information of the progenitor, we calcu¬ 
lated and plotted the lower limit of the bolometric luminos¬ 
ity as the function of the effective temperature and plotted 
them on Fig. 1131 It shows that assuming that the progenior 
was an RSG star, we should have seen it in filter F^lAW, 
if its initial zero-age main sequence mass exceeded 14 Mq 
and both in hlters F606W and FSldlF if its mass was 
Mzams > 20Mq. Note, however, stellar evolution track 
calculations have several uncertainties. The final luminosity 
of a star with a given mass depends heavily on the treatment 
of the mixing processes, e.g. a low er overshooting parame - 
ter leads to lower luminosity. As lEldridee fc ToutI (|2004| ~1 
discussed, by turning the convective overshooting off in the 
STARS code, the final luminosity of a star with a certain ini¬ 
tial mass increases. Therefore the upper limit of the progen¬ 
itor mass in the case of SN 2009ib w ould be 2 — 3 higher 
witho ut overshooting. In addition, IWalmswell fc Eldrideel 
(l2012li discussed that during their lives red supergiant stars 
(which were found to be the progenitors of several SNe II-P) 
produce dust that is destroyed in the SN explosion, there¬ 
fore the extinction before the explosion can be significantly 
higher than the one estimated after the explosion, also lead¬ 
ing to higher upper mass limits. Both of these uncertainties 
point toward a higher mass limit, therefore it is reasonable 
to consider that the upper mass limit of an RSG star too 
faint to be undetected in the pre-explosion images would be 
between 14 — 17 Mq 

In summary, there are several possible interpretations 
of the pre-explosion HST images. These scenarios imply that 



100 hff 

Figure 13. A HR diagram showing the evolution tracks of mas¬ 
sive stars with solar-like metallicity from the STARS code. The 
filled circle marks the position of the yellow source “A”. It lies 
on the evolution track of an Mzams = 20 Mq star. The empty 
circle marks the position of an RSG progenitor which was calcu¬ 
lated based on the assumption that only the FSIAW band de¬ 
tection belongs to the progenitor (see text). This star lies on the 
Mzams = 16 Mq track. The bolometric luminosities calculated 
from the measured magnitude upper limits are plotted as a func¬ 
tion of the effective temperature with black dashed {FA50W), 
red (F606VF) and blue (F814VF) continuous lines. This figure 
shows, that if the progenitor were an RSG star more massive 
than ~ 14 Mq, we would have been able to detect it with the 
filter F814VF and if more than ~ 20 Mq, both with FdOdW and 
F814VF. 


the progenitor of SN 2009ib had an initial mass between 14 — 
20 Mq. However, high resolution, late time HST images are 
necessary to be able to r ule out some of the se possibilities, 
similarly as was done bv iMaund et al.l (l2015l ). 
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6.2 Hydrodynamical modelling 

We constrain the main physical properties of SN 2009ib at 
the explosion (namely the ejected mass, the progenitor ra¬ 
dius and the explosion energy) throngh the hydrodynamical 
modeling of the SN observables (i.e. bolometric light curve, 
evolution of line velocities and continuum temperature at 
the photosphere). 

We use the same approach adopted for other SNe 
(e.g. 2007od, 2009E, 2009N, 2012A, and 2012a w; see 


Inserra_et_^ ]|201]JjPastoreno_et^jj201^ jTakats et al.ll2014l : 


Tomasella et al. 120131 : DaH’Ora et al.l^l4l . respectivelvl. in 

which a simultaneous fit of the observables against model 
calculations is performed and two codes are employed for 
computing the models. The first one is the semi-a nalytic 
code described in detail by IZampieri et al.l ll2003l l which 
is used to perform a preparatory study aimed at individ¬ 
uating the approximate location in the parameter space 
describing the SN progenitor at the explosion. The re¬ 
sults of such study are utilised to guide the more realis¬ 
tic, but time consuming model calculations performed with 
the second code which is the general-relativisti c, radiation- 
hydro dyna mics Lagrangian code pr esented in IPumo et al.l 
ll2010ll and Purno_&_^am£ier (l2oill l . Its main features are 
(cfr. also IPumo fc Zampieo 2013l l: a) a fully implicit La¬ 
grangian approach to the solution of the coupled non-linear 
finite difference system of relativistic radiation-hydro equa¬ 
tions, b) an accurate treatment of radiative transfer in all 
regimes (from the one in which the ejected material is opti¬ 
cally thick up to when it is completely transparent), and c) a 
description of the evolution of ejected material which takes 
into account both the gravitational effects of the compact 
remnant and the heating effects linked to the decays of the 
radioactive isotopes synthesized during the SN explosion. 

Based on the adopted explosion epoch (MJD= 
55041.3; Sect. ld.lTTTl . bolometric luminosity and nickel mass 
(M(®®Ni) = 0.046 Mq; Sect. [5]), the best fit gives values 
of total (kinetic plus thermal) energy of 0.55 foe, progeni¬ 
tor radius of 2.8 x lO*^® cm (~ 400 Rq) and envelope mass 
of 15 Mq (see Fig. I14II . The estimated uncertainty on the 
modelling parameters are about 15 per cent. 

Adding the mass of the compact remnant (~ 1.5 — 
2.0 Mq) to that of the ejected material, we estimate the 
mass of the progenitor of SN 2009ib at the explosion as 
16.5 — 17 Mq, which can be even higher if we consider the 
mass that the progenitor may have lost during its life. This 
estimate is consistent with the mass range of RSG pre¬ 
cursors, and somewhat higher than the upper limit of the 
progenitor mass determined from the pre-explosion images 
(Sect. [6TT| . This is in agreement with results in the litera¬ 
ture that found that progenitor mass estimate from hydro¬ 
dynamical modelling are in general la r ger than those from 
direct detections (see e.g. ISmar ti l2009l : iMaguire et al.|[201^ : 
iTomasella et al.ll2013l ). 


6.3 Comparison to non-LTE models 

We compared the observational properties of SN 2009ib 
to the no n-LTE time-dependen t radiative-transfer simu- 
lations of Dessart et al] ll2013ll . They used MESA STAR 
dPaxton et al.l 20 111 . 2013l l to evolve a 15 Mq main-sequence 
star. Changing parameters such as mixing length, overshoot, 



phase [day] 



phase [day] 


Figure 14. Comparison of the evolution of the main observables 
of SN 2009ib with the best-fit models computed with the general- 
relativistic, radiation-hydrodynamics code (total energy 0.55 foe, 
pre-explosion radius 2.8 X 10^® cm, envelope mass 15 Mq). Top, 
middle, and bottom panels show the bolometric light curve, the 
photospheric velocity, and the photospheric temperature as a 
function of time. To better estimate the photosphere velocity from 
observations, we use the minima of the profile of the Sen lines 
which are considered good tracer of the photosphere velocity in 
Type II SNe. For the sake of completeness, the best-fit model 
computed with the semi-analytic code (total energy 0.4 foe, ini¬ 
tial radius 2.4 X 10^® cm, envelope mass 12.4 Mq) is also shown. 


rotation and metallicity, they created a grid of pre-SN model 
stars. Then they generated a p i ston- driven explosion with 
vId llLivnelfl99d: Dessart et al.||2010| ) and evolved th e SN 


with CMFGEN llHillier k. Miller 1998: ^ssart fc Hillierll20^ : 


IPessart et al.l l2008l : iHillier fc Dessart 201^ un til late into 
its nebular phase (for details see lPessart e t al.ll20 13l). The 
model properties can be found in Table 1 of Dessart et al.l 
1I2OI3I) . They compared the light curve and spectra of 
SN 1999em to the models, and found the only model able 
to reproduce the colour evolution is the one called ml5mlt3, 
the one where the RSG progenitor have significantly smaller 
pre-explosion radius (~ 500 Rq) than the rest. 

We compared the properties of SN 2009ib to two of 
the models: to ml5mlt3 mentioned above, and to the model 
ml5e0p6, which had lower ejecta kinetic energy than the 
baseline model (0.6 foe instead of 1.2 foe). Fig. [TJ] com¬ 
pares the light curve of SN 2009ib to the two models, while 
Fig. ll6l shows the spectral comparison. It is clear from these 
comparisons, that while the model ml5mlt3 fits better the 
colour evolution and the length of the plateau phase, the 
model ml5e0p6 reproduces better the brightness and the 
spectral evolution of the SN. Consequently, the best-fitting 
model might be the combination of these two, a model with 
a smaller progenitor radius and small ejecta kinetic energy. 
The progenitor radius of the ml5mlt3 model and the kinetic 
energy of the ml5e0p6 model are close to the values we ob¬ 
tained from the hydrodynamical modelling. The hydrogen 
envelo pe mass for both models is 10.2 Mq. IPessart et al.l 
1I2OI3I) notes, that the ejecta mass depends also on the He 
core mass, which only can be estimated via nebular phase 
modelling, which goes beyond the scope of this paper. There 
is a broad agreement between the our early nebular phase 
spectra and the models, therefore the best-fitting model 
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Figure 15. Comparison of the lig ht curve of SN 2009 ib to 
the models ml5mlt3 and ml5e0p6 of lPessart et all ll2013h . The 
ml5e0p6 model is scaled to the distance of the SN and matches 
its brigthness well, while the brightness of the model ml5e0p6 is 
higher, therefore we applied an additional vertical shift of 0.88 
mag for better comparison. 


should not be very far from those we examined here, but 
a more careful comparison to a broader grid of models is 
necessary to draw more accurate conclusions. 


6.4 Comparison to other SNe II-P 


The brightness of SN 2009ib places it among the “interme¬ 
diate” luminosity SNe ITP such as SNe 2008in and 2009N, 
however it also shows many differences. SN 2009ib has an 
unusually long plateau, which can be explained by its rela¬ 
tively massive hydrogen-rich envelope. The hydrodynamical 
modelling revealed an ejecta mass of ~ significantl y 

higher than those of SN 2009N (1 1.5 M 0 lTakhts et, alll2ni4l l. 
SN 2008in 11,8.5 M 0 l^iro et al.lli^'l . 

The estimated amount of ®®Ni mass (0.046 Mq) is closer 
to those of the normal SNe ITP, e.g. 200 4et (0.056 M 0 , 
iMaguire et al.l[2010l ~l or 1999em ('0.036 M 0 . IUtrot^l20oit i. 
and significantly higher than those of t he subluminous 
SNe such as SN 2005cs (0.003 — 0.004 M 0 . IPastorello et H] 
[ 200 ^ and intermediate luminosity SNe (e.g. 0.020 Mq 
for SN 2009N, iTakats et al.|[20l3 ). The explosion energy 
(0.55 foe) on the other hand is similar to those of the in¬ 
termediate luminosity SNe (e.g. 0.48 foe for SN 2009N) and 
between those o f the normal ones (e.g. 1.3 foe for SN 1999em, 
IZamnieril l2007ll . and those of the subluminous ones (e.g. 

O. 2 - 0.4 foe for SN 2005cs). 

In Fig. [ 17 ] we compare the physical parameters of 
SN 2009ib and its progenitor to those of other SNe IT 

P. The majority of the data we used here were pub¬ 
lished in IZampie^ (l2007lj who carried out hydrodynam¬ 
ical modelling to estimate the parameters of SNe and 
compared them using similar plots as Fig. 1171 We also 


inclu ded in t he comparison SN e 2007od ( Inserra^_et_^ 


201 J), 2009N llTakats et al.l l2014h . 2 009bw ( IrLserra^.et_aI 


2012li . 1 995ad, 2009dd, and 2010aj lllnserra et al.l l2013lj. 


2012A (iTomasella et ahl l2013lj. 2012 aw ( DaU’Ora et al.l 
l2014h . 2005cs and 2008in ( Sniro et al.ll20l3 l. The physical 
parameters of these objects were estimated with the same 
hydrodynamical modelling technique that we use in this pa¬ 


per (Sect. 16.211 . In the Figure we highlighted the sublumi¬ 
nous SN 2005cs, the normal SN 1999em and the peculiar 
SN 1987A, which objects are often used to represent their 
subgroups, as well as SNe 2008in and 2009N, which have the 
same brightness as SN 2009ib. 

Fig. ll7l con6rms the intermediate nature of SN 2009ib, 
which falls mostly in the middle of the examined parameter 
space. Its nickel mass is among the higher ones in this sam¬ 
ple. Its envelope mass is the closest to that of the peculiar 
SN 1987A, which emerged from a blue supergiant progenitor 
and has different observational properties than the majority 
of SNe II-P. The luminosity as well as the explosion energy 
of SN 2009ib is very close to those of SNe 2008in and 2009N, 
but the nickel mass and hydrogen envelope mass are higher 
than those of the other two, which can explain the longer 
plateau phase. 


7 SUMMARY 

In this paper we present data and analysis of SN 2009ib in 
NGC 1559. This object lies between the luminosity ranges of 
subluminous and normal SNe ITP, belonging to the growing 
sample of the so-called “intermediate” luminosity objects, 
such as SNe 2008in and 2009N. However, this object shows 
significant differences from those, further blurring the line 
between the sub-groups of SNe ITP. 

Among the most important observational properies of 
SN 2009ib are its moderate bolometric luminosity, unusually 
long plateau (lasting 120 days after explosion) and mod¬ 
erate expansion velocities. The spectra is very similar to 
that of the subluminous SN 2002gd, but shows differences 
from those of other subluminous and intermediate luminos¬ 
ity SNe. Comparing the NIR spectra to those of other ob¬ 
jects, we find that the Sr ll lines weaker and the C l lines are 
significantly stronger than usual (noting, however, that the 
available sample is quite limited). 

We estimated the distance to SN 2009ib using both the 
expanding photosphere method and the standardized candle 
method. We also determined the distance to SN 1986L, an¬ 
other SN ITP that exploded in the same galaxy. Combining 
these results with distances from the literature we obtained 
the value of U = 19.8 ± 3.0 Mpc (/r = 31.48 ± 0.31 mag). 
Using this distance, we constructed the bolometric light 
curve of the SN, and from the tail luminosity we esti¬ 
mated the mass of the ®®Ni produced in the explosion as 
0.046 ± 0.015 Mq, a value comparable to those of normal 
SNe ITP, such as SN 2004et. 

We examined the pre-explosion images of the galaxy 
taken by the HST, searching for the progenitor. We found 
two sources close to the SN, one of them coinciding with 
its position. We detected this source in the filters F606W 
and TSldVr. Assuming it is a single star, the yellow colour 
and the luminosity of the source would indicate a type Cl 
star with Mzams ~ 20 Mq. It is also possible, however, that 
it is a blend of two or more stars, and the yellow colour is 
due to the superposition of the light from those. Assum¬ 
ing that in the F&Q&W band an accompanying blue star is 
detected while the F’814VT band detection belongs to the 
progenior and this progenitor is an RSG star, we obtained 
the mass estimate Mzams ~ 16 Mq. In addition, we exam¬ 
ined the possibility that the progenitor of the SN is too faint 
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Figure 16. Comparison of the optical spectra of SN 2009ib with those of the models ml5mlt3 (left) and ml5e0p6 (right) of lPessart et alj 
1120131 1. There is an obvious colour mismatch in the case of the day 16 spectrum with both models. The overall fit is not bad, though the 
velocities and the strength of the features do not match everywhere. 


to appear in HST images. Assuming that it exploded as a 
RSG star, we found that the upper limit of this mass to be 
14- 17Mq. 

Using hydrodynamical modelling, we estimated the pro¬ 
genitor parameters at the explosion, and found the ejecta 
mass 15 Mq, the total explosion energy 0.55 foe and the 
inital radius as 2.8 x 10^® cm. Adding the mass of a compact 
remnant, the estimated progenitor mass is 16.5 — 17 Mq. 

We compared the light curve and spectra of SN 2009ib 
to the non-LTE time-depen dent radiative-transfer simula¬ 
tions of dPessart et ahll^Old ). where the main sequence mass 
of the model star was 15 Mq . We found that the models with 
reduced inital radius match reasonably the colour evolution, 
plateau length and spectral evolution of the SN, while we 
need the model with smaller kinetic energy to match the 
absolute brightness. 

In summary, we found that SN 2009ib may have 
emerged from a red supergiant star with a main sequence 
mass of 14 — 20 Mq and an initial radius of ~ 400 — 500 Mq . 
Its explosion had moderate energy of ~ 0.6 foe, resulting in 
a SN with intermediate luminosity and expansion velocity, 
but somewhat elevated amount of ®®Ni and an unusually 
long plateau phase. 


ACKNOWLEDGEMENTS 

We thank the referee for the useful comments that helped to 
improve this manuscript. We thanks Luc Dessart for sending 
us his models in digital format and for his comments; Mark 
Phillips for sending us the photometric data of SN 1986L; 
and Giorgos Leloudas for observing SN 2009ib. 

K.T. was supported by the Gemini-CONIGYT Fund, 
allocated to the project N° 32110024 and by CONIGYT 
through the FONDEGYT grant 3150473. Support for K.T., 
G.P., F.B. and M.H. is provided by the Ministry of Economy, 
Development, and Tourism’s Millennium Science Initiative 
through grant IC12009, awarded to the Millennium Insti¬ 
tute of Astrophysics, MAS. We acknowledge the TriGrid VL 
project and the INAF-Astronomical Observatory of Padua 
for the use of computer facilities. M.L.P. acknowledges fi¬ 
nancial support from GSFNSM and from the PRIN-INAF 
2011 “Transient Universe: from ESO Large to PESSTO” 
(P.I. S. Benetti). S.B., N.E.R. and E.C. are also partially 
supported by the same PRIN-INAF. N.E.R. acknowledges 
the support from the European Union Seventh Frame¬ 
work Programme (FP7/2007-2013) under grant agreement 
n. 267251 ’’Astronomy Fellowships in Italy” (AstroFIt). M. 
F. received support from the European Union FP7 pro¬ 
gramme through ERG grant number 320360. S.J.S. acknowl- 


© RAS, MNRAS 000,[T]-?? 

























SN 2009ib in NGC 1559 17 






Figure 17. Comparison of the explosion energy (E), nickel mass (Mni) and ejected envelope mass (Menv) of SN 2009 ib (red diamond ) 
estimated via hydrodynamical modelling to those of other SN e IhR_The_^ources of_^e dataMrsed_for cornparison_areJZamDier|j[2003), 
Inserr^^taJjj[201^ , llnserra et al.l ||2013 1 , llnserra et al.l 1I2OI3I I , iTomasella et all ll2013h , iTakats et al.l ll2014l l , iDall’Ora et alJ ll2014l ~) md 


Spiro^ all i2014h 7 In the bottom left figure the size of the points are inversely proportional to log(i?), to better show their position 


along that axis. 
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APPENDIX A: PHOTOMETRIC TABLES OF 
SN 2009IB 

In this Section we include the the magnitudes of the local 
sequence stars that we used to calibrate the photometric 
measurements of SN 2009ib fTables [All IA21IA3II . and the 
calibrated light curves of the SN fTables [A4IIA51IA6II . All 
these tables are available in a machine-readable format in 
the online version of the paper. 
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Table Al. UBVRI magnitudes of the local sequence of stars used for the calibration. 


Star 

«2000 

<52000 

U 

B 

V 

R 


I 

1 

04:17:07.209 

-62:46:13.44 


19.367 (0.078) 

17.850 (0.038) 

16.830 (0.031) 

15.690 

(0.024) 

2 

04:17:05.345 

-62:45:39.15 


16.663 (0.038) 

15.754 (0.030) 

15.212 (0.022) 

14.757 

(0.011) 

3 

04:17:28.126 

-62:45:47.45 


18.610 (0.042) 

17.569 (0.027) 

16.997 (0.044) 

16.526 

(0.023) 

4 

04:17:38.715 

-62:45:32.07 


18.608 (0.040) 

17.116 (0.044) 

16.203 (0.045) 

15.359 

(0.016) 

5 

04:17:04.664 

-62:45:52.59 


17.047 (0.042) 

16.437 (0.027) 

16.094 (0.018) 

15.733 

(0.011) 

6 

04:18:02.142 

-62:45:30.44 


17.632 (0.034) 

16.641 (0.038) 

16.066 (0.024) 

15.542 

(0.014) 

7 

04:18:01.888 

-62:45:58.22 


15.674 (0.017) 

14.989 (0.026) 

14.601 (0.024) 

14.234 

(0.018) 

8 

04:18:06.207 

-62:45:57.24 


18.252 (0.049) 

17.759 (0.033) 

17.489 (0.042) 

17.151 

(0.050) 

9 

04:18:02.508 

-62:47:20.28 


17.508 (0.032) 

16.218 (0.028) 

15.412 (0.034) 

14.758 

(0.010) 

10 

04:18:11.929 

-62:48:49.89 


16.690 (0.025) 

15.569 (0.034) 

14.831 (0.033) 

14.258 

(0.022) 

11 

04:18:06.386 

-62:49:30.78 


19.033 (0.056) 

18.045 (0.029) 

17.354 (0.029) 

16.811 

(0.027) 

12 

04:17:51.943 

-62:50:34.43 


13.703 (0.017) 

13.038 (0.019) 

12.671 (0.025) 

12.332 

(0.014) 

13 

04:17:40.451 

-62:50:46.79 


13.818 (0.020) 

13.300 (0.018) 

12.996 (0.014) 

12.705 

(0.012) 

14 

04:17:29.532 

-62:50:08.59 


15.561 (0.014) 

14.948 (0.018) 

14.590 (0.022) 

14.249 

(0.016) 

15 

04:17:32.022 

-62:49:05.93 


18.418 (0.037) 

17.592 (0.022) 

17.100 (0.026) 

16.679 

(0.023) 

16 

04:17:06.264 

-62:48:48.01 


16.909 (0.013) 

16.161 (0.033) 

15.692 (0.020) 

15.266 

(0.014) 

17 

04:17:07.593 

-62:47:35.37 


18.212 (0.028) 

17.386 (0.037) 

16.891 (0.035) 

16.434 

(0.018) 

18 

04:17:05.610 

-62:47:15.30 


18.809 (0.023) 

18.216 (0.060) 

17.819 (0.020) 

17.471 

(0.037) 

19 

04:17:53.330 

-62:46:48.26 

16.524 (0.101) 

15.666 (0.015) 

14.713 (0.008) 

14.462 (0.238) 

13.753 

(0.026) 

20 

04:17:44.525 

-62:45:42.38 

19.372 (0.192) 

19.199 (0.085) 

18.466 (0.054) 

18.127 (0.019) 

17.655 

(0.040) 

21 

04:17:46.119 

-62:47:44.06 

20.591 (0.459) 

19.173 (0.012) 

18.105 (0.016) 

17.502 (0.018) 

16.860 

(0.020) 

22 

04:17:40.856 

-62:48:20.15 

20.120 (0.047) 

19.709 (0.050) 

18.977 (0.054) 

18.579 (0.052) 

18.087 

(0.081) 


Table A2. u'g'r'i'z' 

magnitudes of the local sequence of stars used for the calibration. 





Star 

0^2000 

<52000 

u’ 

1 


I 


i 

i’ 

2 


1 

04:17:07.209 

-62:46:13.44 



18.525 

(0.099) 

17.284 

(0.034) 

16.310 

(0.020) 

15.865 

(0.064) 

2 

04:17:05.345 

-62:45:39.15 



16.216 

(0.020) 

15.468 

(0.023) 

15.215 

(0.021) 

15.097 

(0.035) 

3 

04:17:28.126 

-62:45:47.45 



18.104 

(0.047) 

17.296 

(0.032) 

17.002 

(0.026) 

16.878 

(0.046) 

4 

04:17:38.715 

-62:45:32.07 



17.831 

(0.072) 

16.604 

(0.034) 

15.942 

(0.033) 

15.611 

(0.059) 

5 

04:17:04.664 

-62:45:52.59 



16.729 

(0.021) 

16.308 

(0.023) 

16.163 

(0.018) 

16.129 

(0.027) 

6 

04:18:02.142 

-62:45:30.44 



17.142 

(0.021) 

16.338 

(0.034) 

16.019 

(0.032) 

15.868 

(0.046) 

7 

04:18:01.888 

-62:45:58.22 



16.673 

(0.057) 

15.310 

(0.019) 

14.843 

(0.026) 

14.678 

(0.026) 

8 

04:18:06.207 

-62:45:57.24 



18.008 

(0.035) 

17.701 

(0.028) 

17.574 

(0.016) 

17.590 

(0.052) 

9 

04:18:02.508 

-62:47:20.28 



16.859 

(0.036) 

15.743 

(0.031) 

15.279 

(0.034) 

15.050 

(0.046) 

10 

04:18:11.929 

-62:48:49.89 



16.145 

(0.019) 

15.140 

(0.031) 

14.740 

(0.027) 

14.557 

(0.036) 

11 

04:18:06.386 

-62:49:30.78 



18.561 

(0.042) 

17.648 

(0.038) 

17.282 

(0.034) 

17.133 

(0.041) 

12 

04:17:51.943 

-62:50:34.43 

14.689 

(0.043) 

13.326 

(0.031) 

12.906 

(0.016) 

12.761 

(0.020) 

12.730 

(0.019) 

13 

04:17:40.451 

-62:50:46.79 

14.699 

(0.052) 

13.528 

(0.033) 

13.216 

(0.007) 

13.124 

(0.017) 

13.122 

(0.021) 

14 

04:17:29.532 

-62:50:08.59 

16.514 

(0.039) 

15.235 

(0.019) 

14.822 

(0.004) 

14.683 

(0.018) 

14.645 

(0.020) 

15 

04:17:32.022 

-62:49:05.93 



17.958 

(0.062) 

17.376 

(0.017) 

17.134 

(0.020) 

17.065 

(0.043) 

16 

04:17:06.264 

-62:48:48.01 



18.097 

(0.074) 

16.516 

(0.040) 

15.947 

(0.017) 

15.719 

(0.016) 

17 

04:17:07.593 

-62:47:35.37 



17.825 

(0.056) 

17.143 

(0.024) 

16.902 

(0.010) 

16.797 

(0.033) 

18 

04:17:05.610 

-62:47:15.30 



18.508 

(0.052) 

18.070 

(0.032) 

17.925 

(0.037) 

17.905 

(0.054) 
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Table A3. Magnitud es of the local sequence of stars used for the calibration of NIR photometry taken from the 2MASS catalogue 
l|Skrutskie et al^2006h . 


Star 

^2000 

<52000 

J 

H 

A 

04:17:36.46 

-62:42:41.79 

11.723 (0.027) 

11.216 

(0.024) 

B 

04:17:43.44 

-62:43:57.30 

12.409 (0.027) 

12.051 

(0.024) 

C 

04:18:01.92 

-62:45:58.35 

13.703 (0.030) 

13.331 

(0.030) 

D 

04:17:53.33 

-62:46:48.26 

13.068 (0.026) 

12.572 

(0.026) 

E 

04:18:02.56 

-62:47:20.28 

13.945 (0.032) 

13.219 

(0.030) 

F 

04:18:11.95 

-62:48:50.11 

13.399 (0.027) 

12.768 

(0.026) 

G 

04:17:51.92 

-62:50:34.94 

11.826 (0.026) 

11.500 

(0.024) 

H 

04:17:40.48 

-62:50:46.96 

12.303 (0.028) 

12.010 

(0.024) 

I 

04:17:29.58 

-62:50:08.81 

13.790 (0.027) 

13.417 

(0.030) 

J 

04:17:19.90 

-62:47:51.73 

11.848 (0.027) 

11.570 

(0.023) 
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Table A4. UBVRl magnitudes of SN 2009ib. 


Date 

MJD 

U 

B 

V 

R 

I 

Telescope 

2009/08/11 

55054.4 


16.030 

(0.016) 

15.854 (0.015) 

15.612 (0.014) 

15.499 (0.014) 

PROMPT 

2009/08/13 

55057.3 


16.241 

(0.157) 

15.810 (0.263) 

15.617 (0.271) 

15.491 (0.181) 

NTT 

2009/08/20 

55063.3 


16.431 

(0.042) 

15.875 (0.031) 

15.547 (0.033) 

15.409 (0.059) 

CTIO 

2009/08/21 

55064.4 


16.569 

(0.020) 

15.920 (0.015) 

15.575 (0.015) 

15.439 (0.013) 

PROMPT 

2009/08/23 

55066.3 


16.657 

(0.053) 

15.931 (0.065) 

15.562 (0.049) 

15.381 (0.061) 

CTIO 

2009/08/23 

55066.4 


16.694 

(0.020) 

15.964 (0.015) 

15.599 (0.013) 

15.420 (0.013) 

CTIO 

2009/08/24 

55067.4 


16.748 

(0.027) 

15.975 (0.015) 

15.584 (0.013) 

15.410 (0.013) 

PROMPT 

2009/08/25 

55068.4 


16.780 

(0.025) 

15.986 (0.014) 

15.632 (0.015) 

15.418 (0.015) 

PROMPT 

2009/08/26 

55069.3 


16.771 

(0.023) 

15.942 (0.042) 

15.555 (0.047) 

15.373 (0.056) 

CTIO 

2009/08/27 

55070.4 


16.863 

(0.029) 

16.019 (0.015) 

15.649 (0.015) 

15.415 (0.013) 

PROMPT 

2009/08/29 

55072.3 


16.957 

(0.026) 

16.041 (0.017) 

15.675 (0.014) 

15.430 (0.015) 

PROMPT 

2009/08/29 

55072.4 


16.891 

(0.043) 

15.975 (0.029) 

15.592 (0.047) 

15.377 (0.024) 

CTIO 

2009/08/31 

55074.3 


17.014 

(0.020) 

16.076 (0.016) 

15.649 (0.013) 

15.426 (0.013 

PROMPT 

2009/09/01 

55075.4 


16.985 

(0.047) 

16.036 (0.030) 

15.596 (0.022) 

15.371 (0.015) 

CTIO 

2009/09/03 

55077.3 


17.090 

(0.034) 

16.110 (0.018) 

15.685 (0.015) 

15.401 (0.022) 

PROMPT 

2009/09/05 

55079.4 


17.146 

(0.035) 

16.123 (0.025) 

15.664 (0.013) 

15.403 (0.038) 

PROMPT 

2009/09/08 

55082.3 


17.147 

(0.099) 

16.116 (0.054) 

15.636 (0.051) 

15.377 (0.039) 

CTIO 

2009/09/09 

55083.3 


17.218 

(0.023) 

16.121 (0.014) 

15.696 (0.013) 

15.396 (0.013) 

PROMPT 

2009/09/10 

55084.3 


17.215 

(0.031) 




PROMPT 

2009/09/11 

55085.3 


17.259 

(0.034) 

16.122 (0.015) 

15.684 (0.014) 

15.391 (0.013) 

PROMPT 

2009/09/11 

55085.4 


17.286 

(0.041) 

16.134 (0.024) 

15.639 (0.025) 

15.324 (0.015) 

CTIO 

2009/09/12 

55086.3 


17.293 

(0.018) 

16.121 (0.014) 

15.675 (0.013) 

15.373 (0.013) 

PROMPT 

2009/09/14 

55088.4 


17.301 

(0.024) 

16.126 (0.025) 

15.690 (0.019) 

15.356 (0.031) 

PROMPT 

2009/09/14 

55088.4 


17.248 

(0.042) 

16.086 (0.029) 

15.632 (0.033) 

15.362 (0.022) 

CTIO 

2009/09/16 

55090.3 


17.316 

(0.024) 

16.157 (0.016) 

15.720 (0.015) 

15.372 (0.013) 

PROMPT 

2009/09/18 

55092.3 


17.330 

(0.045) 

16.178 (0.033) 

15.758 (0.071) 

15.364 (0.102) 

PROMPT 

2009/09/20 

55094.3 




16.184 (0.016) 

15.726 (0.015) 

15.388 (0.014) 

PROMPT 

2009/09/22 

55096.3 


17.402 

(0.022) 

16.221 (0.015) 

15.727 (0.015) 

15.398 (0.013) 

PROMPT 

2009/09/22 

55096.3 


17.389 

(0.046) 

16.179 (0.053) 

15.658 (0.028) 

15.385 (0.033) 

CTIO 

2009/09/24 

55098.3 


17.423 

(0.034) 

16.219 (0.014) 

15.720 (0.013) 

15.396 (0.012) 

PROMPT 

2009/09/26 

55100.3 


17.438 

(0.031) 

16.216 (0.018) 

15.739 (0.016) 

15.384 (0.015) 

PROMPT 

2009/09/30 

55104.3 


17.499 

(0.061) 

16.195 (0.024) 

15.708 (0.021) 

15.332 (0.020) 

CTIO 

2009/10/05 

55109.3 


17.516 

(0.022) 




PROMPT 

2009/10/07 

55111.3 


17.545 

(0.021) 




PROMPT 

2009/10/09 

55113.3 


17.506 

(0.027) 

16.241 (0.015) 

15.730 (0.013) 

15.347 (0.013) 

PROMPT 

2009/10/12 

55116.3 


17.523 

(0.038) 

16.216 (0.016) 

15.728 (0.013) 

15.362 (0.012) 

PROMPT 

2009/10/13 

55117.3 


17.528 

(0.109) 




PROMPT 

2009/10/17 

55121.3 




16.204 (0.014) 

15.703 (0.013) 

15.344 (0.012) 

PROMPT 

2009/10/19 

55124.4 

19.425 (0.349) 

17.531 

(0.157) 

16.213 (0.145) 

15.715 (0.243) 

15.323 (0.106) 

NTT 

2009/10/27 

55131.3 


17.506 

(0.040) 

16.206 (0.020) 

15.683 (0.017) 

15.328 (0.016) 

PROMPT 

2009/10/30 

55134.3 




16.204 (0.026) 



PROMPT 

2009/11/03 

55138.2 


17.601 

(0.020) 

16.255 (0.015) 

15.731 (0.013) 

15.344 (0.012) 

PROMPT 

2009/11/07 

55142.3 


17.576 

(0.030) 




PROMPT 

2009/11/08 

55143.3 




16.192 (0.018) 

15.711 (0.015) 

15.359 (0.013) 

PROMPT 

2009/11/12 

55147.2 


17.583 

(0.039) 




PROMPT 

2009/11/13 

55148.2 




16.262 (0.017) 

15.729 (0.015) 

15.384 (0.015) 

PROMPT 

2009/11/16 

55151.1 


17.645 

(0.029) 

16.302 (0.014) 

15.758 (0.013) 

15.399 (0.012) 

PROMPT 

2009/11/22 

55157.2 

19.741 (0.455) 

17.749 

(0.188) 

16.369 (0.184) 

15.874 (0.201) 

15.457 (0.179) 

NTT 

2009/11/23 

55158.1 





15.812 (0.019) 

15.448 (0.018) 

PROMPT 

2009/11/24 

55159.2 




16.365 (0.020) 

15.806 (0.013) 

15.458 (0.012) 

PROMPT 

2009/12/02 

55167.1 




16.478 (0.017) 



PROMPT 

2009/12/03 

55168.1 


17.891 

(0.039) 




PROMPT 

2009/12/07 

55172.1 


17.950 

(0.062) 

16.602 (0.018) 

15.983 (0.017) 

15.632 (0.014) 

PROMPT 

2009/12/09 

55174.1 


18.054 

(0.050) 




PROMPT 

2009/12/10 

55175.2 


18.039 

(0.080) 

16.697 (0.018) 

16.071 (0.013) 

15.721 (0.013) 

PROMPT 
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Table A4 - continued 


Date 

MJD U B 

V 

R 

I 

Telescope 

2009/12/12 

55177.1 

18.152 (0.064) 




PROMPT 

2009/12/14 

55179.2 

18.376 (0.061) 




PROMPT 

2009/12/15 

55180.1 

18.393 (0.059) 




PROMPT 

2009/12/16 

55181.1 


17.098 (0.024) 

16.394 (0.017) 

15.999 (0.015) 

PROMPT 

2009/12/18 

55183.0 

18.692 (0.065) 

17.298 (0.030) 



PROMPT 

2009/12/22 

55187.3 


17.570 (0.022) 

16.809 (0.021) 

16.402 (0.015) 

PROMPT 

2009/12/24 

55189.0 


17.680 (0.023) 

16.864 (0.027) 


PROMPT 

2009/12/24 

55190.2 

19.266 (0.260) 

17.768 (0.137) 

16.986 (0.250) 

16.559 (0.159) 

NTT 

2010/01/04 

55200.1 

19.206 (0.109) 




PROMPT 

2010/01/07 

55203.0 


17.931 (0.090) 



PROMPT 

2010/01/08 

55204.0 



17.035 (0.036) 


PROMPT 

2010 /01/12 

55208.1 

19.250 (0.168) 




PROMPT 

2010/01/17 

55213.1 

19.282 (0.121) 




PROMPT 

2010/01/18 

55214.2 


18.003 (0.032) 



PROMPT 

2010/01/24 

55221.2 


17.999 (0.190) 

17.140 (0.284) 

16.636 (0.199) 

NTT 

2010/04/17 

55304.0 


18.702 (0.221) 

18.143 (0.373) 

18.207 (0.231) 

NTT 


Telescopes: PROMPT - Panchromatic Robotic Optical Monitoring and Polarimetry Telescopes; CTIO - CTIO 1.3 m SMARTS 
telescope + ANDICAM; NTT - 3.6-m New Technology Telescope + EFOSC2. 
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Table A5. u'g'r'i'z' magnitudes of SN 2009ib obtained with the PROMPT telescopes. 


Date 

MJD 

u' 

9' 

r' 

i' 



2009/08/11 

55054.4 

16.148 (0.035) 

15.921 (0.015) 

15.781 (0.013) 

15.869 (0.015) 

15.935 

(0.019) 

2009/08/21 

55064.4 

17.651 (0.069) 

16.222 (0.015) 

15.797 (0.013) 

15.840 (0.015) 

15.836 

(0.018) 

2009/08/23 

55066.4 

17.791 (0.084) 

16.293 (0.015) 

15.828 (0.012) 

15.840 (0.014) 

15.820 

(0.017) 

2009/08/24 

55067.4 

18.022 (0.074) 

16.314 (0.015) 

15.830 (0.012) 

15.838 (0.026) 



2009/08/25 

55068.4 


16.352 (0.017) 





2009/08/26 

55069.4 

18.204 (0.087) 






2009/08/27 

55070.4 


16.420 (0.019) 





2009/08/28 

55071.3 

18.306 (0.080) 






2009/08/29 

55072.3 


16.458 (0.016) 





2009/08/30 

55073.4 

18.517 (0.068) 






2009/08/31 

55074.3 


16.486 (0.015) 





2009/09/01 

55075.3 

18.600 (0.088) 






2009/09/05 

55079.4 


16.548 (0.019) 

15.910 (0.015) 

15.850 (0.018) 

15.788 

(0.053) 

2009/09/09 

55083.4 


16.582 (0.014) 

15.912 (0.013) 

15.854 (0.013) 

15.765 

(0.015) 

2009/09/11 

55085.3 


16.594 (0.014) 

15.913 (0.013) 

15.858 (0.014) 

15.770 

(0.015) 

2009/09/13 

55087.4 


16.635 (0.014) 

15.910 (0.013) 

15.847 (0.015) 

15.766 

(0.059) 

2009/09/14 

55088.4 





15.757 

(0.026) 

2009/09/15 

55089.3 


16.654 (0.016) 





2009/09/16 

55090.4 





15.750 

(0.021) 

2009/09/17 

55091.3 


16.663 (0.014) 

15.924 (0.013) 

15.864 (0.014) 



2009/09/18 

55092.4 




15.831 (0.026) 



2009/09/19 

55093.2 


16.693 (0.015) 

15.921 (0.012) 

15.854 (0.012) 

15.747 

(0.014) 

2009/09/21 

55095.3 


16.720 (0.017) 

15.940 (0.016) 

15.853 (0.021) 

15.752 

(0.018) 

2009/09/23 

55097.4 


16.721 (0.017) 

15.943 (0.013) 

15.852 (0.014) 

15.753 

(0.015) 

2009/09/25 

55099.3 


16.735 (0.016) 

15.953 (0.013) 

15.866 (0.013) 

15.745 

(0.015) 

2009/10/01 

55105.3 


16.766 (0.014) 

15.955 (0.016) 

15.867 (0.018) 

15.740 

(0.029) 

2009/10/05 

55109.3 


16.774 (0.016) 





2009/10/07 

55111.3 


16.784 (0.015) 





2009/10/08 

55112.3 


16.783 (0.016) 

15.965 (0.013) 

15.865 (0.015) 

15.725 

(0.017) 

2009/10/12 

55116.3 


16.798 (0.016) 

15.975 (0.015) 

15.864 (0.013) 

15.696 

(0.015) 

2009/10/18 

55122.3 



15.972 (0.011) 

15.855 (0.012) 

15.695 

(0.014) 

2009/10/24 

55128.2 


16.804 (0.014) 





2009/10/28 

55132.3 


16.807 (0.014) 

15.976 (0.015) 

15.847 (0.012) 

15.695 

(0.015) 

2009/11/02 

55137.2 


16.811 (0.014) 

15.974 (0.013) 

15.841 (0.013) 

15.704 

(0.017) 

2009/11/05 

55140.3 


16.826 (0.015) 

15.979 (0.015) 

15.843 (0.017) 

15.696 

(0.025) 

2009/11/08 

55143.2 


16.825 (0.018) 





2009/11/11 

55146.2 


16.840 (0.018) 





2009/11/13 

55148.2 


16.865 (0.016) 





2009/11/16 

55151.2 



16.006 (0.013) 

15.881 (0.019) 

15.740 

(0.019) 

2009/11/24 

55159.2 



16.051 (0.015) 

15.940 (0.016) 

15.761 

(0.015) 

2009/12/03 

55168.1 


17.140 (0.020) 





2009/12/07 

55172.2 


17.259 (0.017) 

16.229 (0.013) 

16.146 (0.013) 

15.927 

(0.019) 

2009/12/09 

55174.2 


17.326 (0.019) 





2009/12/10 

55175.2 


17.371 (0.026) 





2009/12/11 

55176.2 




16.428 (0.032) 



2009/12/12 

55177.2 


17.474 (0.022) 





2009/12/14 

55179.2 


17.603 (0.021) 





2009/12/15 

55180.2 


17.689 (0.021) 





2009/12/18 

55183.1 


17.985 (0.030) 





2009/12/20 

55185.1 


18.198 (0.030) 





2009/12/21 

55185.3 



16.981 (0.024) 

16.854 (0.028) 

16.556 

(0.032) 

2009/12/23 

55188.1 


18.398 (0.030) 





2009/12/24 

55189.1 


18.395 (0.030) 





2009/12/26 

55191.1 



17.184 (0.022) 

17.039 (0.028) 

16.715 

(0.024) 

2009/12/27 

55192.1 



17.222 (0.021) 

17.077 (0.028) 

16.735 

(0.028) 

2009/12/28 

55193.2 



17.316 (0.025) 

17.119 (0.030) 

16.768 

(0.035) 

2010 /01/01 

55197.1 


18.419 (0.025) 





2010 /01/02 

55198.1 



17.361 (0.021) 

17.176 (0.020) 

16.878 

(0.020) 

2010/01/04 

55200.1 


18.494 (0.024) 





2010/01/13 

55209.0 


18.498 (0.032) 





2010/01/18 

55214.1 


18.576 (0.035) 





2010/02/13 

55240.1 


18.742 (0.027) 
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Table A6. JH magnitudes of SN 2009ib. 


Date 

MJD 

J 

H 

2009/08/12 

55056.3 

15.36 (0.22) 

15.08 

(0.56) 

2009/08/20 

55063.3 

15.15 (0.27) 

15.05 

(0.67) 

2009/08/25 

55068.3 

15.13 (0.18) 

14.85 

(0.63) 

2009/08/30 

55073.4 

15.06 (0.11) 

14.66 

(0.24) 

2009/09/10 

55084.2 

14.91 (0.15) 

14.57 

(0.45) 

2009/09/15 

55089.3 

14.92 (0.07) 

14.51 

(0.56) 

2009/09/21 

55095.2 

14.87 (0.16) 

14.48 

(0.12) 

2009/09/26 

55100.2 

14.84 (0.22) 

14.51 

(0.64) 

2009/10/12 

55116.3 

14.67 (0.68) 

14.45 

(0.11) 

2009/10/17 

55121.3 

14.77 (0.18) 

14.45 

(0.26) 

2009/10/23 

55128.2 

14.73 (0.05) 

14.43 

(0.35) 

2009/11/06 

55141.2 

14.76 (0.17) 

14.47 

(0.37) 

2009/11/11 

55146.3 

14.72 (0.12) 

14.50 

(0.36) 
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